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Abstract 

This letter is concerned with transmit and receive filter optimization for the K-user MIMO inter- 
ference channel. Specifically, linear transmit and receive filter sets are designed which maximize the 
weighted sum rate while allowing each transmitter to utilize only the local channel state information. 
Our approach is based on extending the existing method of minimizing the weighted mean squared 
error (MSE) for the MIMO broadcast channel to the K-user interference channel at hand. For the case 
of the individual transmitter power constraint, however, a straightforward generalization of the existing 
method does not reveal a viable solution. It is in fact shown that there exists no closed-form solution for 
the transmit filter but simple one-dimensional parameter search yields the desired solution. Compared 
to the direct filter optimization using gradient-based search, our solution requires considerably less 
computational complexity and a smaller amount of feedback resources while achieving essentially the 
same level of weighted sum rate. A modified filter design is also presented which provides desired 
robustness in the presence of channel uncertainty. 
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I. Introduction 

To achieve high spectral efficiency, much effort has been focused on improving the achievable 
rate of multiple-input multiple-output (MIMO) interference channels [l]-[3]. A notable scheme 
in this area, the interference alignment (IA) technique of [4] confines all undesired interferences 
from other communication links into a pre-define subspace and achieves a maximum-capacity 
scaling. However, it is also known that IA can only offer a suboptimal sum rate at finite signal- 
to-noise ratios (SNRs) [3]. 

In this letter, we aim at maximizing the sum rate in the K-user MIMO interference channel. We 
consider two linear transceiver design methods. One is for the sum-power-usage-limit constraint 
and the other applies to the per-transmit-node power-usage constraint. The former can be viewed 
as a network-level constraint whereas the latter is more of a device-level constraint. In both 
designs, to maximize the weighted sum rate (WSR), we pursue minimization of the weighted 
mean squared error (WMSE). The idea of maximizing the WSR via receiver-side WMSE 
minimization was originally developed for the multi-user MIMO broadcast channel [5]. Our 
sum-power-constrained method could be seen as a generalization of the approach of [5] to cover 
the K-user MIMO interference channel and can be obtained as a direct extension of the method 
in [5]. However, our individual-power-constrained method is not a direct generalization of the 
method of [5] due to multiple power constraints. In fact, unlike in the case of the broadcast 
channel, we show that there is no closed-form solution for the minimum WMSE transmit filter, 
although a simple one-dimensional search for the power-adjusting parameter leads to the desired 
solution. Using simulation results and analysis, we verify that both proposed schemes achieve 
the maximum WSR with lower computational complexity than the gradient-based optimization 
of the transmit and receive filters [2]. Also, unlike in [2], [4], [6], our schemes require only the 
local channel state information (CSI) (i.e., each transmitter needs to know only the CSI of the 
links originating from itself whereas the MIMO interference channel precoder designs in [2], 
[4], [6] require the CSI for all links). Additionally, we discuss modified transceiver design that 
provides significant robustness in the presence of inaccurate CSI. 

Related ideas for the MIMO interference channel can also be found in [3], [6]-[10]. In [3], 
[8], the minimum MSE (MMSE) transceiver is designed without considering different weights 
for the MSEs at multiple receivers. In [6] suboptimal MSE weights are used. In contrast, our 
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weighted MMSE transceiver design relies on a set of MSE weights that provides a direct link 
between the weighted MMSE (WMMSE) and WSR criteria. The WMMSE-based weighted utility 
maximization is also considered in [7], but there only a single data stream is assumed between 
a given user pair. A very similar idea on maximizing WSR via WMSE minimization under the 
individual power constraint has been discussed in [9]. But, unlike in our approach, the inter- 
dependency between the transmit-power-adjusting Lagrange multiplier and the precoding matrix 
has not been considered in [9]. In our individual-power-constrained transceiver design, this inter- 
dependency is handled by introducing one-dimensional search for the Lagrange multiplier. This 
means that the method of [9] requires recursive optimization based on exchanges of filter- setting 
information among all transmitters. Our method does not require recursive filter adjustment and 
no data exchanges are needed among transmitters. 1 Finally, we present a modified transceiver 
design method for the imperfect-CSI environment and analyze the computational complexity as 
well as the required feedback amount in comparison with the gradient descent method [2]. 

The following notations are used. We employ upper case boldface letters for matrices and 
lower case boldface for vectors. For any general matrix, X, X T , X*, ~X H , Tr(X), det(X), 
vec(X), SVD(X) denote the transpose, the conjugate, the Hermitian transpose, the trace, the 
determinant, the stack columns, and the singular value decomposition of X, respectively. The 
symbol || ■ || 2 indicates the 2-norm of a vector. The symbol I n denotes an identity matrix of size 
n. 

II. System Model 

We consider the MIMO interference channel where precoding can only be done over one 
transmission slot. As shown in Fig. 1, K source nodes simultaneously transmit independent 
data streams to their desired destination nodes and generate co-channel interference to all other 
undesired nodes. In this system each source node {S^} is equipped with M antennas and each 
destination node {D fc } has iV antennas (k e {1 ~ K})- The MIMO channels from Sj to 
Dj are modelled by H,; e C NxM G {1 ~ K}) whose coefficients are independent and 

'The independently conducted and recently published work of [10], which was brought to our attention by an anonymous 
reviewer, also pursues maximization of the WSR via weighted MSE minimization. The transceivers in [10] do become the same 
as our proposed individual-power-constrained transceivers when each base station serves a single user. Relative to the work of 
[10], this letter includes the sum-power-constrained method as well as a method to handle mismatched CSI. 
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identically distributed (i.i.d) complex Gaussian random variables with CAf(0, erf). We assume 
that the channel information is only locally available, i.e., each node knows only the coefficients 
for the channel link originating from itself. Note that the precoder designs of [2], [4], [6] are 
based on the availability of the global channel information. Let e C dxl denote the symbol 
vector from S k with E[s fe s^] = l d where d is the number of data streams for D fc , d < M,N 
and the value of d is chosen to meet the feasibility of degree of freedom [11]. Also \ k e C Mxd 
denotes the precoding matrix for S k . Then, the N x 1 received signal vector at D k is represented 
as 

K 

y k = HfefcVfeSfe + ^ HfciVjSj + n fc , (1) 

where n k denotes the i.i.d complex Gaussian noise vector at D k with CAf(0, ct^Ijv)- Then, D k 
combines its received signal with XJ k £ C dxN to decode the desired signals: 

K 

Sfc = U fe y fc = U k H kk V k s k + U k ^2 HfejVjSi + U fc n fe . (2) 

i^k 

Our goal is to find {V fe } and {U k } that maximize the WSR under the sum-power constraint and 
also the individual-power constraint. We assume a unit noise variance (o 2 n = 1) without losing 
generality. 

III. Weighted Sum Rate Maximization 
First consider finding {V fc } that maximizes 

K 

J^^Rk subject to ^ Tr(V fe Vf) = P T or Tr(V fe Vf ) = P k Mk (3) 

k=l k 

where the subscript k points the source node and its intended destination node, fj, k denotes 
the weight, R k is the achievable rate, P T represents the maximum sum power allowed for all 
transmitters and P k is the A;-th node's maximum transmit power. With Gaussian signaling, the 
achievable rate takes the well-known form: 

R k = log {det(Tv + S^HfcfcVfcVf H&) }, (4) 

where = 1 N + Yli^k HfciVjV^H^. We attempt to solve this WSR maximization problem 
by minimizing the weighted receiver MSE, as has been done for the MIMO broadcast channel 
[5]. This approach was also attempted for the K-user MIMO interference channel in [9] under 
the individual-power constraint, but our solution is different as elaborated below. 
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A. Relationship between achievable rate and error covariance matrix 

To understand the link between the WSR maximization problem and the WMSE minimization 
problem in the K-user MIMO interference channel, we need to clarify the relationship between 
the achievable rate and the error covariance matrix. This argument is parallel to one given in [5] 
for the MIMO broadcast channel. For the MMSE receive filter at D fc , we write 

^(MMSE) =argmin ]E||XJ A ,y A , — S fc j || 

K 

=VfHf fc (^H fci V,VfH^ + I N )-\ (5) 

i=i 

and the error matrix for D fe is given by 

E fc =E{ ( ur £) y, - s.xur^v* - s*n 

=(I w + $ fc 1 H fcfe V fc VfH^)- 1 . (6) 

Comparing (4) and (6), the relationship between the achievable rate and the error covariance 
matrix is established as: 

R k = \og {de^E, 1 )} (7) 

which, not surprisingly, is identical to the relationship between the rate and the error covariance 
matrix for the case of the MIMO broadcast channel [5]. Apparently, though, the error covariance 
matrix E fc here is different from that of the broadcast channel due to the presence of multiple 
sources. Note that this relationship between the achievable rate and the error covariance matrix 
holds for any {V fc }, implying that (7) is true with either transmit power constraint. 

B. MSE weight design 

Now consider finding {V fc } that solves the following WMMSE problem: 

K 

min^Tr(W fc E fc ) subject to ^ Tr(V fc Vf) = P T or Tr(V fe Vf ) = P k Vfc, (8) 

k=l k 

where W fc e C dxd represents the MSE weight. Again following the argument of [5], the MSE 
weights can be chosen so that both WSR and WMMSE problems have a common solution. For 
this, set up the Lagrangians for (3) and (8): 

K K K 

Cwsr = ~Y, ^Rk + 9\(J2 Tr(V fe Vf ) -F r ) + (l-9)(^2 A fe (Tr(V fe Vf ) - P fe )) 

k=l k=l k=l 
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and 

K K K 

Cwmse = ^Tr(W fc E fe ) + #A(^Tr(V fe Vf ) - P T ) + (1 - 6) ( ^ A fc (Tr(V fe Vf ) - P fc )) 
fc=i fc=i fc=i 

respectively, where 6> selects the desired power constraint ('# = V for the sum power constraint 
and '9 = 0' for the individual power constraint), A and {X k } denote the Lagrange multipliers for 
the two transmit power constraints. Next, equate their gradients obtained via the matrix derivative 
formulas: rf{ln(det(X))} = Tr{X-M(X)}, d{Tr(X)} = Tr{d(X)}, vec{d(X)} = d{vec(X)}, 
Tr(X T Y) = vec(X) T vec(Y). Subsequently, the resulting MSE weight can be found as 

W k = -^Er 1 . (9) 
k ln(2) k V ' 

Note that the choice of the MSE weights {W fe } is irrelevant to the transmit power constraint, 

which makes sense as {W fc } are receiver-side design parameters. 

C. Sum power constrained precoder design 

We are now ready to find the transmit precoding matrix that minimizes the WMSE under the 
sum-power constraint, i.e., find {V fe } that minimizes 

K 

^E[Tr{W fe (s fc - /T^kXs* - r%) H }] subject to ^Tr(V fc Vf ) = P T (10) 

k=l k 

where {W fc } is set according to (9) and f3 is a scaling parameter. With matrix derivative formulas, 
the WMMSE transmit filter that satisfies (10) can be shown to be 

v fc = /?v;, (id 

where V k = (* k + S^^f h M ) "H&UfWfc, * fe = Eii HgUf W,U,H ife , and (3 = 



£ fc Tr(V fc V fe ) 

This result is a rather straightforward generalization of the WMMSE precoder in the broadcast 
channel. It can indeed be seen that setting H ki = H kk for all i, our solutions (5), (9), and (11) 
reduce to the respective receive filter, MSE weight and transmit filter solutions obtained for the 
multi-user MIMO broadcast channels through WMSE minimization [5]. 
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D. Individual power constrained transceiver design 

Now let us consider the individual-power-constrained network. We proceed to find the transmit 
filter that minimizes the weighted MSE: 

K 

^E[Tr{W fc (s fc -s fe )(s fc -s fc )"}] subject to Tr(V fe Vf) = P k Vfc. (12) 
k=\ 

Again equating the gradients of the Lagrangians corresponding to the WMMSE and WSR 
maximization procedures and using the matrix derivative formulas, the WMMSE transmit filter 
at S fc is found as: 

V, = + A fe I M ) ~ V fc Uf W fc (13) 

where X k is set to satisfy the transmit power constraint at S k and again {W fc } are as given in 
(9). Unlike the sum-power-constrained WMMSE precoders of (11), for which the power control 
parameters are found in closed form, here we resort to a numerical method to find X k , due to 
the inter-dependency between V k and X k in (13). Fortunately, based on the following lemma, 
\ k can be found with simple one-dimensional (1-D) numerical search. 

Lemma 1: The per-node transmit power, Tr(V fc (Afc)V fe (Afc)^), is a monotonically decreasing 
function of X k . 

Proof: Let SVD(\I> fc ) = Q k Ti k Q, k . Then, the transmit power at S k is given by 

Tr{(V fc (A fc )V fc (A fc ) H } =Tr{(Q fc S fe Qf + A J M )- 1 Hf fe Uf W fe U fe H fcfe (Q fe S fc Qf + A^m)" 1 } 

=Tr{(S fe + A fc I M r 2 QMuf W k U k H kk Q k } 

= E 7 [ ^fy 2 , 04) 

{<Tk,i + Afc) 2 

where II fc = Q k ll kk XJ k W k \J k 'H kk Q, k , [Ilfe]^ is the ii-th element of and a k>i is the i-th 
element of S fc . Because X k > 0, ^{V^X^V^Xk) 11 } is monotonically decreasing with X k . ■ 
Note that the proper set of MSE weights for the K-user MIMO interference channel has 
already been derived in [9] in the process of establishing a connection between the WMMSE 
problem and the WSR maximization problem. In [9], though, the transmitter S k is expressed as 
a function of itself as well as transmitters at the other nodes, i.e. V k = f({Vi, • • ■ , Vjf}). The 
consequence of this formulation is that the transmitter solution in [9] cannot be found without 
recursive calculation and additional filter- setting information exchanges among all transmit nodes. 
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In contrast, our transmit filter design is based on a clear recognition of the inter-dependency 
between X k and V fc , and as a result the proposed transmit filter (13) can be found through a 
simple 1-D numerical search with no additional information ({V^} (I ^ k)) exchanges needed 
among the transmit nodes. 

E. Iterative algorithm to maximize the weighted sum rate 

In the previous sections, we found the MSE weights and then subsequently WMMSE receive 
and transmit filters with both the sum-power constraint and the individual-power constraint. Each 
of three sets of parameters - MSE weights, transmit filters and receive filters - is derived assuming 
the other sets are given. In practice, to find optimum WSR solutions, the inter-dependencies 
between the parameters are handled with the following iterative or alternating optimization 
algorithm. 

Algorithm 1 Obtaining the optimal WSR transceivers via the WMMSE criterion 
Initialize I = and {Vj^}, calculate Flfum- 
repeat 

l:=l + l 

Step 1: Calculate Ujj. |{vf _1) } for all k using (5). 
Step 2: Calculate wf|{vf _1) } for all k using (9). 

Step 3: Calculate V^|{uf } }, {Wj°} for all k using (11) for the sum power constrained 
case or (13) for the individual power constrained case, 
until \Rsum — Rsum\ < e, where e is some arbitrarily small value and R sum = J2 k UkRk- 



The algorithm is common to both the sum-power-constrained design and the individual-power- 
constrained design. This algorithm is provably convergent to a local optimum; this can be shown 
by proving monotonic convergence of an equivalent optimization problem based on expanding the 
WSR maximization problem of (3) to add the MMSE weights and receive filters as optimization 
variables, as has been done for the MIMO broadcast channel in [5]. We note, however, that this 
algorithm does not guarantee the global optimal solution, since the WMMSE minimization (8) is 
not jointly convex over all input variables. To reasonably approach the optimal solution one must 
resort to repeated runs of the algorithm using different initial settings, or, for computationally 
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efficient initialization, choose {Vj^} in Step 1 from the right singular matrices of {H fefe } or 
from random matrices generated according to the normal distribution with zero mean and unit 
variance [8]. 

IV. Robust transceiver design for imperfect channel information 

In practical scenarios, mismatch between the true channel {H^} and the estimated channel 
(denoted by {H^}) is inevitable because of the channel estimation errors [12]. In this section, 
we design robust transceivers for mitigating the performance degradation caused by channel 
mismatch. We assume that {H^} is related to {H^} by Hjj = Hj,- + Ay where the elements 
of are independent and identically distributed (i.i.d.) complex Gaussian random variables 
with variance a\ [12]. Then, the received signal can be rewritten as 

K 

s k = \J k {H kk - A kk )V k s k + U fe ^(Hk - Afe^VjSj + \J k n k (15) 

where {V k } and {U k } are computed from {H^} with no knowledge of the presence of {A^}. 
We try to mitigate the effect of channel mismatch by minimizing the appropriate metrics averaged 
over Aj/s. 

1 ) Modified MSE weight: Following the design procedure in previous sections, a modi- 
fied version of the MMSE receiver filter is found as \J k = Vf Hf fc (ES= i H^Vf + 
E*=i <riTr(Av.)Iiv + Ijv) -1 , where SVD(V;Vf) = QjA^Qf . The modified MSE weights 
that force the optimum solutions of the WSR maximization and WMMSE problems to be 



identical are derived as W k = I ^J|yE fc 1 , where E fc = (I N + <£ fc ^^VfeVf H^)^ 1 and $ fe = 

Iiv + Eik HwViVf Hjg + Ef=i ^(A^)l N . 

2 ) Robust transceiver design with the sum power constraint: The modified transmit filters are 
derived based on the following optimization problem: 

K 

min^E[Tr{W fc (s fe - p-%)(s k - p-%) H }] subject to ^Tr(V fc Vf ) = P T . (16) 

k=l k 

Utilizing matrix derivative formulas, the resultant modified-WMMSE transmit filters are obtained 

as 

V fc = pV k , (17) 

where V; = (^+ ^^ C '^ ) Im+E^i ^Tt(A 0< )I m ) ~H*Uf W fc , = y/=^ y 
*k = E^HlUfW.U.H^, and SVD(UfW i U i ) = Q^. 
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3 ) Robust transceiver design with the individual power constraint: The optimization problem 
to derive the modified precoder is 

K 

min^E[Tr{W fe (s fe -s fe )(s fe -s fe )"}] subject to Tr(V fc Vf) = P k Wk. (18) 

k=l 

With the matrix derivative formulas, the modified-WMMSE transmit precoder at S k with the 
individual power constraint is written as 

V k =(^ k + \ k l M + J2^A^ii^M) Hf fe UfW fc (19) 

i=i 

where the power control parameter \ k is also found by numerical 1-D search. 

Note that, for the above derivations, we have assumed that the value of the channel error 
variance a\ is perfectly known. In the practical systems, the channel error variance can be 
estimated through an appropriate statistical approach [13]. Below, we also present numerical 
performance results corresponding to the cases where the error variance is not perfectly known. 

V. Discussion: Computational complexity, channel state information 

In this section, we analyze computational complexity and required feedback resources. For 
comparison, we also analyze those of the gradient descent method of [2]. 

A. Computational complexity 

We consider the number of complex multiplications as a complexity measure. As summarized 
in the Table I, the number of complex multiplications is proportional to the number of itera- 
tions. The proposed method with the sum-power constraint which has a single iteration loop 
is computationally the most efficient. Whereas both the proposed method with the individual- 
power constraint and the gradient descent method require double iteration loops, i.e., the outer 
loop for updating the sum rate and the inner loop for adjusting the Lagrange multiplier (in 
the case of the proposed method) or for updating the step size (in the case of the gradient- 
based method). Calculating the gradient and adjusting the step size require more computational 
resources. According to simulation, when SNR = 10 dB which is in the mid SNR regime, K = 4, 
M = N = 5, and d = 2, the minimum average numbers of iteration for the convergence of sum 
rate, updating the step size of gradient method and 1-D search with bisection method are 10, 10 
and 10, respectively. In accordance with these, ii = 10, J 2 = 10 and I 3 = 10 are chosen. The 
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symbols c]y, c 2 NM and c 3 N denote the computational complexity of a matrix inversion of iV x N 
matrix, a singular value decomposition ofNxM matrix, and a Cholesky factorization of N x N 
matrix, respectively. The corresponding values for those variables are f A^ 3 , 7NM 2 + AM 3 , and 
|iV 3 , respectively [14]. Fig. 2 shows comparison when M = N = 5 and d = 2 2 . As expected, 
for the same WSR values the proposed method with the sum-power constraint has the least 
complexity while the gradient descent algorithm is the most computationally complex. 

B. The amount of required feedback information 

To find the optimized transmit precoders, each transmit node requires feedback information. 
As illustrated in Table II, feedback information is composed of CSI and coefficients for filter 
updating. For a given transmission slot, CSI feedback is required once, but the filter coefficients 
are updated several times due to the iterative optimization algorithm. Although the proposed 
method requires a larger amount of feedback information for the iteratively updated coefficients 
such as MSE weights {W fe } and receive filter coefficients {Ufc} than the gradient descent method 
does, the amount of CSI feedback for the proposed method is smaller than for the gradient descent 
method. This is because, unlike the global CSI requirement of the gradient-based method, the 
proposed methods need only local CSI. From Table II, we observe that as the network size 
grows (i.e., K increases) the required feedback resources for local CSI and coefficient updating 
increase linearly, but those for global CSI increases quadratically. Fig. 3 clearly shows that 
with Ii = 10 the proposed methods are advantageous in terms of required feedback resources, 
especially for larger K. Note that, for the transmit power adjustment, the sum-power-constrained 
method additionally requires iterative update of the scalar parameter Tr{^._^ fc V-V- H }, but the 
size of this parameter is negligible compared to other matrix parameters. 

VI. Numerical Results 

In this section, we provide the numerical results related to the WSR performances. The SNR for 
the sum-power-constrained network, SNR = and that for the individual power constrained 
network, SNR fc = ^PsVfc, are derived assuming P T = K, P k = lVfe and a 1 — 1, i.e., 

2 To see the effect of the number of K, we fixed M = N = 5, even though the degree of freedom (DoF) is not achievable 
when K > 5 
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SNR = SNR fc = a\. The results are averaged over 1000 independent trials. Fig. 4 shows 
the average WSR performance of the proposed methods for M = N = 5 (when K = 4), 
M = N = 6 (when K = 5), and d = 2. For fairness, all schemes are initialized with the right 
singular matrices of the intended channels. For the sum-power constraint, we set the weights 
to be ii-i — 2 and fi k = 0.25 (k ^ 1), which were chosen rather arbitrarily except that ^ is 
made considerably larger than (j, k to bring out the performance advantage of the sum-power 
constraint. The performance of the sum power constraint method should be better than that of 
the individual power constraint method because the former, which is less stringent, is able to 
allocate more power to the higher weighted transmitter to maximize the WSR. When the weights 
are equal, /j, k — 1 Vfc, the performance of both proposed schemes and that of the conventional 
gradient descent method are nearly identical. Note that, as explained in section IV, the proposed 
methods achieve these performances with less computational complexity and a smaller amount 
of feedback resources than the gradient descent method. Compared to the performance of the 
MMSE transceiver without the MSE weights [3], [8] (curves labelled "Simple MMSE"), the 
advantage of designed MSE weights is clearly shown as SNR grows. Fig. 5 demonstrates the 
effectiveness of the robust design with either transmit power constraint in presence of channel 
uncertainty when K = 4 and a\ = O.lcr^. As SNR grows, the amount of leakage interference 
due to CSI imperfection also increases. This is why the performance is saturated in the high 
SNR regime in Fig. 5. To reflect a potential error in estimating a\, we model the channel 
error variance as o\ + of, where o\ is the actual channel error variance and of indicates over- 
estimation. As shown in Fig. 5, at SNR = 15 dB at most 3 % sum rate losses are shown when 
of = 0.1(7^. Although not shown, same results were observed for under-estimating the channel 
estimation error variance. 

VII. Conclusion 

In this letter, we have studied a linear transceiver design method for the K-user MIMO 
interference channel. To maximize the weighted sum rate with less computational complexity 
and a smaller amount of feedback resources, the proposed transceivers are designed in the 
weighted MMSE sense with suitably chosen MSE weights. Also, the proposed transceiver design 
considers both the sum-power-usage constraint and the individual-power constraint. Through 
numerical simulation, we have demonstrated that the weighed- sum-rate performances of the 
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proposed schemes approach that of the existing gradient descent method. The proposed methods 
have clear advantage in terms of processing requirements as well as feedback resources over the 
gradient-based technique. Also, modified versions of proposed schemes have been provided for 
compensating channel mismatch. 
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Num. of users, K 



Fig. 2. Complexity comparison between the gradient-descent-based method and the proposed methods 
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Num. of users, K 



Fig. 3. The amount of feedback information at each source node to design precoder 



— $ — Grad. descent (no weights) 
Q — Sum pwr. prop. (equal weights) 

- © - Sum pwr. prop. (unequal weights) 
A Ind. pwr. prop. (equal weights) 

-A- Ind. pwr. prop. (unequal weights) 
| Simple MMSE (no weights) 




HC Grad. descent (no weights) 
— — Sum pwr. prop. (equal weights) 
- - Sum pwr. prop. (unequal weights) 
Ind. pwr. prop. (equal weights) 
Ind. pwr. prop. (unequal weights) 
— I — Simple MMSE (no weights) 
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SNR [dB] 



(a) K=4, M=N=5, d=2 (b) K=5, M=N=6, d=2 

Fig. 4. The weighted sum rate performance (equal weight fik = 1, Vfc and unequal weights fii — 2, ji k i — 0.25, k' ^ 1) 
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- © - Ind. pwr. prop. 

Q Ind. pwr. modified (b 2 =0) 
■ -Q- i Ind. pwr. modified ^»0.1o^) 
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— H — Sum pwr. modified (o 2 =0) 
i -X- 1 Sum pwr. modified (o 2 =0.1a 2 ) 
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(b) Unequal weights {p\ = 2, [l% = 0.25, fc 7^ 1) 



Fig. 5. The weighted sum rate performances with imperfect channel state information (<r A = 0.1 x a h , K — 4) 
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(a) Description of each stage for gradient descent method 



STAGE 


Index 


Initialization 


a.l 


Outer loop 


Calculating gradient 


a.2 


Inner loop: calculating step size 


a.3 


Calculating sum rate 


a.4 


Calculating optimal precoders and decoders 


a.5 



(b) Description of each stage for proposed methods 



STAGE 


Index 


Initialization 


b.l 




Calculating the variance of noise and interference 


b.2 




Calculating the receive filter 


b.3 


Loop 


Calculating the error covariance matrix 


b.4 




Calculating the MSE weights 


b.5 




Calculating the transmit filter 
(1-D search is needed for individual power constraint) 


b.6-1 (for sum power constraint) 
b.6-2 (for individual power constraint) 




Calculating sum rate 


b.7 


(c) Number of complex multiplication at each stage 



Index 


Number of complex multiplication 


a.l 


K(M 2 d + 1) + K (K - 1)(1 + 2MNd + N 2 d) + K(2 + 2MNd + N 2 d + N 3 + c^) 


a.2 


h^K(2K - 1)(1 + 2MNd + N 2 d) 
+K(2K - 1)(9 + 2cn + 2MN 2 + 2M 2 N + 2M 2 d + Md 2 )} 


a.3 


h[KI 2 (l2 + l)/2 + KI 2 {2K{K - 1)(1 + 2MNd + N 2 d) 
+2K(2 + 2MNd + N 2 d + N 3 + c^) + K(M 2 d + 1) + 2 + A/d 2 }} 


a.4 


h{K(M 2 d + 1) + K(K - 1)(1 + 2MNd + N 2 d) + K(2 + 2MNd + N 2 d + N :i + <4)} 


a.5 


K(l + 2Md + 2M 2 d + c 2 Md ) + K(K - l)(2MNd + N 2 d) 
+K(2MNd + 2N 2 d + ANd 2 + Md 2 + d :i + + c :i d + c 2 d + c 2 dd ) 


b.l 


K(M 2 d + 1) + K (K - 1)(1 + 2MNd + N 2 d) + K{2 + 2MNd + N 2 d + N 3 + c„) 


b.2 


I\K{K - 1)(1 + 2MNd + N 2 d) 


b.3 


hK(3MNd + 2N 2 d + c„) 


b.4 


hK(2MNd + N 2 d + Nd 2 + c l N + cj) 


b.5 


hKc\ 


b.6-1 


h{K(K - l)(2NMd + Md 2 + M 2 d) + K(Nd 2 + d 3 ) 
+K(3MNd + 2Md 2 + M 2 d + 1 + c l M ) + K(M 2 d + Md) } 


b.6-2 


h {K(K - \){2NMd + Md 2 + M 2 d) + hKM 2 d 
+(7 3 + l)K(3MNd + 2Md 2 + M 2 d + 1 + c^)] 


b.7 


h{K(M 2 d + 1) + K(K -!)(! + 2MNd + N 2 d) + K{2 + 2MNd + N 2 d + N 3 + cjv)} 



TABLE I 

Computational complexity comparison 
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Grad. descent method 


Prop, method 


Global CSI 


Updating coefficients 


Local CSI 


Updating coefficients 


Feedback information 


{H^} 


{V,}, (i + k) 


{H lfc } 


{UJ, {W t } (Ind. pwr.) 


{U,}, , {W,}, J2^ k Tr{vM H } (Sum pwr.) 


Matrix size 


MNK 2 


Md(K - 1) 


MNK 


(Md + d 2 )K (Ind. pwr.) 


(Md + d 2 )K + 1 (Sum. pwr.) 


Feedback resource amount 


MNK 2 + Md(K - l)ii 


MNK + (Md + d 2 )Kh (Ind. pwr.) 


MNK + ((Md + d 2 )K + l)Ji (Sum. pwr.) 



TABLE II 

Summary of required feedback information at the fc-TH transmit node, S^, i, j = 1 ~ K 
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